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ABSTRACT 

We present our analysis of Kepler observations of 29 RR Lyrae stars, based on 138-d of ob- 
servation. We report precise pulsation periods for all stars. Nine of these stars had incorrect 
or unknown periods in the literature. Fourteen of the stars exhibit both amplitude and phase 
Blazhko modulations, with Blazhko periods ranging from 27.7 to more than 200 days. For 
V445 Lyr, a longer secondary variation is also observed in addition to its 53.2-d Blazhko pe- 
riod. The unprecedented precision of the Kepler photometry has led to the discovery of the 
the smallest modulations detected so far. Moreover, additional frequencies beyond the well- 
known harmonics and Blazhko multiplets have been found. These frequencies are located 
around the half-integer multiples of the main pulsation frequency for at least three stars. In 
four stars, these frequencies are close to the first and/or second overtone modes. The ampli- 
tudes of these periodicities seem to vary over the Blazhko cycle. V350Lyr, a non-Blazhko 
star in our sample, is the first example of a double mode RR Lyrae star that pulsates in its 
fundamental and second overtone modes. 
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1 INTRODUCTION 

The ultraprecise photometry by the Kepler space telescope opens 
up the possibility of discovering new phenomen a and shedding 
new l ight on long-standing astrophysical problems JGilliland et al. I 
|201(J) . One of t he most in t erestin g unsolved problems is the physi- 
cal origin of the lBlazhkol I I 19071) effect, an amplitude and/or phase 
modulation of RRLyrae stars. The leading explanatio ns are: (1) an 
obliq ue rotator model that invokes a magnetic field l lShibahashi I 
120000 : (2) a m odel with resonant coupling betwe en radial and non- 
radial mode(s) JDziembowski & Mizerski 120041) ; and (3) a mecha- 
nism invoking a cyclic variati on of the turbule nt convection caused 
by a transient magnetic field dStothers 1120060 . At present none of 
these models explains all the observed properties of stars showing 
the Blazhko effect. It is not even clear whether a modification of the 
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above ideas or new astrophysical processes are needed to solve the 
problem. Comprehensive discussions of the observational a nd the- 
oretica l pr operties of Blaz hko RR Lyrae stars are given by ISzeidll 
dl988h and lKovacs I J2009I) . 

This paper paper describes new properties of RRLyrae stars 
revealed by early data from the Kepler photometer. Here we con- 
centrate on the results based mostly on Fourier analyses. A detailed 
study of all observed stars is beyond the scope of the present paper. 



2 DATA 

A detailed technical d escri ption of the Kepler Mi ssion can be found 
in lKoch et all J2010h and I Jenkins et al. I J2010JlbT) . At the time of 
this writing three long cadence (29.4-min integration time) photo- 
metric data sets have been released to the KASC (Kepler Astero- 
seismic Science Consortium). Altogether 29 RRLyrae stars were 
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Table 1. Basic parameters of 29 RRLyrae stars observed by Kepler (long-cadence) during the runs from QO to Q2. 



KIC 



RA 

[J2000] 



DEC 
[J2000] 



K p 
[mag] 



[d] 



a(P ) Ax 

[10- 5 d] [mag] 



[10- 3 



ID 



Runs 



Note 



magj 



3733346 


19 08 27.23 


+38 48 46.19 


12.684 


0.68204 


1.68 


0.266 


2.2 


3864443 


19 40 06.96 


+38 58 20.35 


15.593 


0.48680 


0.73 


0.331 


2.4 


3866709 


19 42 08.00 


+38 54 42.30 


16.265 


0.47071 


0.85 


0.340 


3.0 


4484128 


19 45 39.02 


+39 30 53.42 


15.363 


0.54787 


1.24 


0.299 


2.9 


5299596 


19 51 17.00 


+40 26 45.20 


15.392 


0.52364 


0.91 


0.195 


1.5 


5559631 


19 52 52.74 


+40 47 35.45 


14.643 


0.62072 


1.48 


0.273 


2.4 


6070714 


19 56 22.91 


+41 20 23.53 


15.370 


0.53410 


0.93 


0.229 


1.7 


6100702 


18 50 37.73 


+41 25 25.72 


13.458 


0.48815 


0.71 


0.206 


1.5 


6183128 


18 52 50.36 


+41 33 49.46 


16.260 


0.56168 


1.07 


0.245 


1.9 


6186029 


18 58 25.69 


+41 35 49.45 


17.401 


0.51293 


1.13 


0.204 


2.0 


6763132 


19 07 48.37 


+42 17 54.67 


13.075 


0.58779 


1.13 


0.280 


2.3 


7030715 


19 23 24.53 


+42 3142.35 


13.452 


0.68362 


1.43 


0.231 


1.8 


7176080 


18 49 24.43 


+42 44 45.56 


17.433 


0.50708 


0.94 


0.357 


3.0 


7198959 


19 25 27.91 


+42 47 03.73 


7.862 


0.56688 


1.26 


0.239 


2.2 


7505345 


18 53 25.90 


+43 09 16.45 


14.080 


0.47370 


1.29 


0.374 


3.4 


7671081 


19 09 36.63 


+43 21 49.97 


16.653 


0.50457 


0.86 


0.313 


2.5 


7742534 


19 10 53.40 


+43 24 54.94 


16.002 


0.45649 


0.77 


0.407 


3.5 


7988343 


19 59 50.67 


+43 42 15.52 


14.494 


0.58115 


1.28 


0.341 


3.0 


8344381 


18 46 08.64 


+44 23 13.99 


16.421 


0.57683 


1.26 


0.322 


2.9 


9001926 


18 52 01.80 


+45 18 31.61 


16.914 


0.55682 


1.13 


0.287 


2.5 


9508655 


18 49 08.37 


+46 1 1 54.96 


15.696 


0.59424 


1.33 


0.339 


3.0 


9578833 


19 09 40.64 


+46 17 18.17 


16.537 


0.52702 


0.99 


0.304 


2.5 


9591503 


19 33 00.91 


+46 14 22.85 


13.293 


0.57139 


1.09 


0.384 


3.2 


9697825 


19 0158.63 


+46 26 45.74 


16.265 


0.55759 


1.06 


0.261 


2.1 


9947026 


19 19 57.96 


+46 53 21.41 


13.300 


0.54859 


0.88 


0.219 


1.6 


10136240 


19 19 45.28 


+47 06 04.43 


15.648 


0.56579 


1.23 


0.270 


2.4 


10789273 


19 14 03.90 


+48 1158.60 


13.770 


0.48029 


0.86 


0.390 


3.4 


11125706 


19 00 58.77 


+48 44 42.29 


11.367 


0.61323 


0.98 


0.179 


1.2 


12155928 


19 18 00.49 


+50 45 17.93 


15.033 


0.43639 


0.71 


0.394 


3.4 



NRLyr 


Q1,Q2 


V2178Cyg 


Q1.Q2 Bl,n 


V715 Cyg 


Q1.Q2 


V8O8 Cyg 


Q1.Q2 Bl 


V782 Cyg 


Q1.Q2 


V783 Cyg 


Q1.Q2 Bl 


V784 Cyg 


Q1.Q2 


ASAS 185038+4125.4 


Q0,Q1,Q2 


V354 Lyr 


Q1.Q2 Bl,n 


V445 Lyr 


Q1,Q2 Bl,n 


NQLyr 


Q0.Q1.Q2 


ASAS 192325+4231.7 


Q0.Q1.Q2 


V349 Lyr 


Q1.Q2 Bl,n 


RRLyr 


Q1.Q2 Bl 


V355 Lyr 


Q2 Bl 


V450 Lyr 


Q1.Q2 Bl 


V368 Lyr 


QLQ2 n 


V1510Cyg 


Q1.Q2 


V346 Lyr 


Q1.Q2 


V353 Lyr 


Q1.Q2 Bl,n 


V350 Lyr 


Q1.Q2 


V366 Lyr 


Q1.Q2 Bl,n 


V894 Cyg 


Q0.Q1.Q2 


V360 Lyr 


Q1.Q2 Bl,n 


V2470 Cyg 


Q0,Q1,Q2 


VI 107 Cyg 


Q1.Q2 n 


V838 Cyg 


Q1.Q2 


ROTSE1 J190058.77+484441.5 


Q0,Q1,Q2 Bl 


VI 104 Cyg 


Q1,Q2 Bl 



observed in this way. One additional R RLyrae star (V355 Lyr) ha s 
been released since the publication by iKolenberg et al. I ( 1201 Oaf) . 
This source was observed as a Director's discretionary targe{j in 
the Kepler Guest Observer Program from quarter 2 onwards, and 
is included here to enrich the KASC sample. Throughout the work 
presented here, we used only the long-cadence datcQ 

The commissioning phase data (Q0) included six of the 
RRLyrae stars between 2009 May 2 to 11 (9.7 d), the observations 
of Ql data began on 2009 May 13 and ended on 2009 June 15 
(33.5 d). The first full quarter of data (Q2) ran from 2009 June 19 
to 2009 September 16 (89 d). At present, the combined number of 
data points for a given star is between 4096 and 6175. Column 8 in 
TableQ]shows the available data sets for each sta£j. 

The telescope is rotated by 90 degrees four times per orbital 
period for best exposure of the solar panels. Accordingly, the Q0 
and Ql data sets were observed at one position and the first roll 
was executed between Ql and Q2. The different apertures applied 
to a given star in the two positions caused the zero points of the raw 
fluxes and the amplitudes of a star to be different for the two rolls. 
We applied simple linear transformations to fit the amplitudes and 
zero points for combining the data. The long time-scale trends were 



1 http://keplergo.arc.nasa.gov/GOprogramDDT.shtml 

2 We have observed several additional RR Lyrae candidates with short ca- 
dence (1-min integration); results from these data will be discussed in a 
forthcoming paper. 

3 Public Kepler data can be downloaded from the web page: 
http://archive.stsci.edu/kepler/ 



removed from the raw fluxe s by a trend filtering algorithm prepared 
for C0R0T RRLyrae data dChadid et al~ll20icin, then fluxes were 
transformed into a magnitude scale, where the averaged magnitude 
of each star was fixed to zero. Measured raw fluxes are in the range 
of 3.1 x 10 10 > F > 1.3 X 10 6 ADU which yields 6 x 10" 6 - 
9 x 10~ 4 mag accuracy for an individual data point. 



3 ANALYSIS AND RESULTS 

As a first step Fourier analyses were performed on the data sets. T o 
this end we used the software packages MuFrAn (Kollath 199i 
Period04 i lLenz & Breger1l2005l) and SigSpec JReegen 11200 
all of which gave similar frequency spectra with similar errors. 
Some basic parameters of the observed stars are summarized in Ta- 
bleQ] The columns of the table shows the ID numbers, J2000 posi- 
tions (RA, DEC), and apparent magnitude {Kepler magnitude K p ), 
all from the Kepler Input Catalog (KIC-lOJj. The next columns 
contain the main pulsation periods and the Fourier amplitude of the 
main frequencies obtained from our SigSpec analysis. These two 
basic parameters (Po and A\) were previously unknown or wrong 
for nine stars (signed by an 'n' in the last column). 

This fact already suggests, that there were hardly any detailed 
investigations on these stars before this work. This is underlined 
by the lack of known stars showing Blazhko effect before Kepler 



4 http://www.konkoly.hu/HAG/Science/index.html 

5 http://archive.stsci.edu/kepler/kiclO/search.php 
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Figure 1. The gallery of Kepler Blazhko stars. The figure shows the complete light curves of eight stars obse r ved with long-cadence d uring the periods Q0 
through Q2. Further light curves are given in Fig.[2]and Fig.[5]and in our parallel papers ISzabo et al. 1120101 ; iKolenberg et al. Il2010bl) . Lines seen running 
through these light curves are visual artefacts caused by beating of the sampling frequency with the pulsation frequencies. They cause no problems in the 
Fourier analysis. *For better visibility the scale of KIC 1 1 125706 is increased by a factor of 1.5. 



Table 2. Period, amplitude and phase properties of the Blazhko stars 



KIC 


GCVS 


Pb 

[d] 


°(Pb) 
[d] 


[mag] 


A0i 


Q 


Addition, freq/ 6 ' 


3864443 


V2178Cyg 


> 200 




> 0.488 


> 0.0014 


0.153 


F2, (PD) 


4484128 


V808 Cyg 


«90 




0.304 


0.012142 


-0.045 


PD 


5559631 


V783 Cyg 


27.7 


0.4 


0.071 


0.001591 


0.156 




6183128 


V354Lyr 


>127 




> 0.245 


> 0.00245 


-0.139 


F2, (F1,PD,F') 


6186029 


V445 Lyr 


53.2( a > 


2.8 


0.968 


0.022442 


0.540 


PD,F1,F2 


7176080 


V349Lyr 


>127 




> 0.060 


> 0.00175 


-0.251 




7198959 


RRLyr 


39.6 


1.8 


0.461 


0.013836 


0.676 


PD 


7505345 


V355 Lyr 


31.4 


0.1 


0.107 


0.003518 


0.136 


PD 


7671081 


V450Lyr 


w 125 




0.391 


0.004882 


0.235 




9001926 


V353 Lyr 


60.0 


7.1 


0.157 


0.004026 


0.033 




9578833 


V366Lyr 


65.6 


2.6 


0.171 


0.003205 


-0.162 




9697825 


V360Lyr 


51.4 


4.3 


0.356 


0.008228 


0.279 


Fl, (PD) 


11125706 




39.4 


2.0 


0.030 


0.001420 


0.540 




12155928 


VI 104 Cyg 


53.1 


0.3 


0.105 


0.002451 


-0.063 





( a 'The star shows a longer time-scale variation than its Blazhko modulation as well. 

'^The pattern of additional frequencies: PD means period doubling; Fl indicates first overtone frequency and its linear combination with fundamental one; 

F2 is as Fl, but with second overtone; F' indicates frequencies with unidentified modes; brackets indicate marginal effects. 



except RRLyr itself. In most of the cases these stars are men- 
tioned in the literature apropos of their discoveries and in some 
cases due to the determination of their position and/or ephemeris. 
The radial velocity measurements of NR Lyr and V894 Cyg were 
used f or kinematic study o f the Galaxy I ILavdenll994l ; lBeers et al. I 
l200rj Ijefferv et al. I l2007h . The Northern Sky Variablity Survey 
(NSVS) included NRLyr, V355Lyr and V2470Cyg. Their NSVS 
light curves - together with mo re than 1 100 other ones - formed the 
basis of the statistical study of iKinemuchi et al. I J2006J) . Th e most 
specifi c inv estigations w ere carried out on V783 Cyg by [Loser] 
( 119790 and ICross I J1991D . who found a period increase between 



1933 and 1990 with the ephemeris of 3D max = 2436394.332 + 
0.62069669-E 1 + 7.5 ■ Kr n .E 2 . According to TableQ]the period is 
still increasing with a good agreement of the rate of Cross's value: 
0.088 ± 0.023 d My" 1. 

The standard errors of the main period and amplitude (in 
columns 6 and 8) were estimated from the accuracy of non-linear 
least-square fits. The last three columns of the table indicate other 
identifications of the stars, the observing runs analysed and the 
presence of a Blazhko effect, respectively. All the periods, ampli- 
tudes and light curve shapes of the 29 stars are typical for RRab 
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Figure 2. Top left: Light curve of V445 Lyr (KIC 6186029). Bottom left: Parts of the light curve from maximum and minimum of the amplitude of the 
Blazhko cycle. Right: Fourier spectrum after the data are prewhitened with the main frequency and its harmonics. The insert is a zoom around the positions of 
the main frequency /n and its first harmonics 2/rj (showed by green arrows) after the 20 highest amplitude frequencies were removed. The complex structure 
of additional frequencies is clearly seen. 



stars pulsating in their radial fundamental mode, therefore, these 
classifications are omitted from Table[T] 



3.1 RR Lyrae stars with amplitude modulation 

Generally, it is an easy task to distinguish amplitude modulated and 
non-modulated Kepler RR Lyrae light curves. A gallery of modu- 
lated light curves is shown in Fig.[T] It is obvious at first glance 
that the modulation cycles are predominantly long and the ampli- 
tude of the effect is clearly visible. Non-sinusoidal envelopes of the 
light curves (see, e.g., V2178 Cyg; KIC 3864443) or moving bumps 
(e.g., V366Lyr; KIC 9578833) are also conspicuous. 

The interesting light curve of V445Lyr (KIC 6186029) is 
shown separately in Fig.[2j The two observed Blazhko cycles are 
surprisingly different. The high amplitude of the Blazhko mod- 
ulation extremely distorts the shape of the light curve. This is 
demonstrated in the small panels of Fig.|2] Signs of complex vari- 
ations are detectable from the Fourier spectrum as well. The spec- 
trum of the data prewhitened with the main pulsation frequency 
and its harmonics shows four peaks around each of the harmon- 
ics (Fig.[3). Two outer peaks at the harmonics can be identified as 
elements of the Blazhko triplets (/n ± /b). Two other side peaks 
closer to the harmonic frequencies show the possible variation of 
the Blazhko effect on a time-scale longer than the observation 
run. This may be a result of a cyclic variation (existence of more 
than one Blazhko modulation), a secular trend, or random changes. 
Several papers have reported mu ltiperiodic and/or unstable be- 
haviours in the Blazhko effect (e.g ., LaCluyze 2004; Collinge et al. 
2006; Kolenberg et al. 2006; Nagv & Kovacs 2006; Sodor et al. 



20061; I Szczygiel & Fabrvckv 1120071 ; IWils et alj|2008l ; ljurcsik et all 



2009c). After investigating the two Blazhko cycles noted here, we 
are not in the position to decide on the nature of this variation, but 
the 3.5-5 year long time base of Kepler's observations will provide 
an excellent opportunity to study this strange behaviour. 

We were systematically searching for low amplitude Blazhko 
RR Lyrae stars. Instrumental trends of the observed fluxes that are 
not properly removed could cause apparent amplitude changes in 
the non-linear magnitude scale. A decreasing trend of the averaged 
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Figure 3. Fourier amplitude spectrum of V445 Lyr (KIC 6186029) around 
the main pulsation frequency /n after the data are prewhitened with the 
main frequency. Beyond elements of the Blazhko triplet /rj ± /b two addi- 
tional peaks are seen (arrows). 



fluxes results in increasing amplitudes in magnitudes and vice- 
versa. Therefore, we always checked the amplitude variation us- 
ing the raw fluxes. We divided the data sets into small sections 
(typically 2-3 d in length), then calculated the amplitude difference 
SAi(t) of the first Fourier component and its averaged value AA\ 
over the whole time span for all sections by a non-linear fit. These 
calculated functions reflect well the variation of pulsation ampli- 
tude seen in the light curves. 

With the help of this tool we found in KIC 11125706 the 
lowest amplitude modulation ever detected in an RR Lyrae star. 
Full amplitude of the maximum light variation A(K p ) max = 
0.015 mag, and the amplitude of the highest side peak in the Fourier 
spectrum is only Ak p (/o + /b) = 0.0022 mag. The lowest pub- 
lished amplit ude of a Blazhko effe ct previously found was the case 
of DM Cyg djurcsik et al.ll2009bl) where A(V) ma * = 0.07 mag, 
Av(/o + /b) = 0.0096 mag and A/ (/ + /b) = 0.0061 mag. The 
two measurements are not strictly comparable because Kepler pass- 
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band is broad in white light. However, th e maximum of its spectral 
response function (see lKoch et al. luOlOf) is about 6000 A between 
Johnson-Cousins filters V and R. 

All the Blazhko RR Lyrae stars found in our sample are listed 
in Table [2] The estimated lengths of the Blazhko cycles are indi- 
cated in the third column. For the shorter periods than the total 
time span they were calculated from the averaged frequency differ- 
ences of the highest side peaks (/ + /b and / — /b), otherwise a 
minimum period is given. The fourth column shows the amplitude 
modulation parameter AAi defined above. 

We note here, the triplet structure has always appeared for all 
Blazhko stars, i.e. no stars show frequency doublets. The amplitude 
of the high frequency peak is higher than the lower for 9 stars. The 
remaining 5 show the opposite pa ttern. The asymmetry parameter 
Q defined bv lAlcock et all J2003h as Q = [A(f + f B ) - A(/ - 
/b)P(/o + /b) + A(f - fs)]' 1 varies from -0.251 to 0.676 
(see the 6th column in Table[2}, however, these values are rather 
preliminary due to the long Blazhko cycles. 

In the past few years, thanks to high precision ground- 
and space-based observations, the known occurrence rate of the 
Blazhko effect among RR Lyrae stars has increased from the for- 
mer estimate of 15-30% to a ratio close to 50% (se e ljurcsik et al. I 
l2009cMChadid et al. l2009llKolenberg et al. |2010J) . It is even pos- 
sible that all RR Lyrae stars may show a Blazhko modulation with 
an increasing frequency of Blazhko stars at lower modulation am- 
plitude. The Kepler measurements provide an ideal tool to test this 
hypothesis. From our AAi values in Table [2] it can be seen that 
we found only two stars with modulation amplitude lower than 
0.1 mag. 

To test our detection limit of the amplitude modulation, 
artificial light curves were generated. Two sets of grids were 
constructed: one for V368Lyr (KIC 7742534) and another for 
KIC 7030715. These stars have the shortest and the longest pul- 
sation periods (0.45649 d and 0.68204 d) in the sample, respec- 
tively. In both cases the Fourier parameters of the main pulsation 
frequency and its significant harmonics were used to build the ar- 
tificial light curves. These were modulated by a simple sinusoidal 
function with amplitudes ranging between 0.1 and 0.001 mag, and 
with modulation periods from 25 to 150 d with a ste p size of 25 d, 
accord ing to the general modulation formula (eq. 2) in lBenko et al. I 
(2009!). Measured averaged fluxes of the non-Blazhko stars are in 
the range of 1.8 x 10 8 > F > 2.7 x 10 6 ADU which means a 
noise between 8 x 10 -5 and 6 x 10 -4 mag. This was taken into 
account by adding Gaussian noise with a — 10 -4 and 5 x 10~ 4 
to the artificial data. The light curves were always calculated at the 
observed points of time. 

In our tests, we reckoned the amplitude modulation as de- 
tectable if the highest Fourier side peak connected to the modu- 
lation exceeds the spectral s i gnific ance (a s ) 5. (For the definition 
of the <t s we refer iReegen I ( 120071) : the correspond ence between 
more popular amplitude signal-to-noise ratio S/N JBreger et al. I 
1 19931) and a s is yielded as a s = 5 « S/N = 3.83). The ob- 
tained limiting values are A(fo + /b) > 0.001-0.002 mag (or 
AAi > 0.005-0.01 mag) depending on the brightness, but highly 
independent of the periods (Po and Pb). Higher sampling rate (i.e. 
short cadence) does not decrease our detection limit because the 
present sampling frequency 48.98 d~ is much higher than a typ- 
ical Blazhko frequency (0.1-0.01 d _1 ), hence each Blazhko cycle 
is covered sufficiently. 

Notwithstanding our efforts, we did not detect any modula- 
tion for 15 RR Lyrae stars in our Kepler sample, however, some 



60 
50 
40 
30 


- 


















20 
10 


















30 
20 


- 


























10 


















































.01 .02 .03 .04 .05 .06 .07 

A(f„+f B ) 

Figure 4. Modulation amplitude distribution of the Blazhko variables in 
the Kepler sample using 0.025 mag (top) and 0.005 mag (bottom) size of 
bins. The modulation amplitude corresponds to the Fourier amplitude of the 
largest amplitude modulation frequency component A(fo + /b)- 



very small amplitude modulation with long period might remain 

undetected. 

Us ing the same number (14) of Blazhko stars ljurcsik et al. I 
d2009a) found an exponential-like distribution of their modulation 
amplitude strengths. Our sample shows a similar behaviour (top in 
Fig.0 when we divide it up in t o the s ame 0.025 mag size of bins 
as were used by ljurcsik et al. I ([2009c). However, the distribution 
seems be to more uniform, when we split our sample into smaller 
size of bins (bottom in Hig.|4jl. Although, our sample is affected by 
small number statistics in this respect, we checked the uniformity of 
the modulation amplitude distribution. We carried out a one-sample 
Kolmogorov-Smirnov (K-S) test over the intervals of (0, AA™ ax ) 
and (0, A(f + /b)" 1 " 1 ), respectively. (The extremely modulated 
star V445 Lyr was omitted from the sample.) Both tests allowed the 
uniform distribution hypothesis with 99 per-cents of probabilities. 
We note that the K-S test uses data directly, without any binning. 



3.2 Phase modulation 

Without any "a priori" knowledge about the nature of the mod- 
ulation, frequency modulation (period changes over the Blazhko 
cycle) and phase modulation can not be distinguished: a detected 
phase variation indicates period changes, and vice- versa. From now 
on we refer to this phenomenon as phase modulation. 

We searched for phase modulation in all of our RR Lyrae stars 
in the same way as was described in the case of amplitude mod- 
ulation. We calculated the Aifii values from the phase variation 
function 8ipi(t). We expressed the phase differences relative to the 
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Figure 5. Left top: Light curves of two Kepler RRLyrae stars, V354Lyr (KIC 6183128) and V360Lyr (KIC 9697825). Right: Fourier spectra after the data 
are prewhitened with the main frequency and its harmonics. Arrows point to the highest peaks connected to additional frequencies. The inserts are a zoom 
around the position of the highest additional peaks. Left bottom: Amplitude of the highest additional peaks vs. time. The typical error is about 0.001 mag - 
smaller than the size of the data points. 



total cycle, that is A0i=A<^i/27r (= SPo/Po). The results can be 
seen in the 5th column of Table[2] 

We detected clear phase modulation for all the studied 
Blazhko RRLyrae stars. The hardest task was to find it in the case 
of V2178 Cyg, where the Blazhko cycle is longer than the data set 
and the phase variation during the observed time span was only 
0.0014 (~1 min). The reality of this small phase variation was 
checked and confirm ed by the sensitive analytical function method 
JKollath et al. 1 120021) . We can detect period variations smaller 
than 1.5 min for three further stars: V783Cyg (KIC 5559631), 
V349Cyg (KIC 7176080), and KIC 1125706. The other extremity 
is represented by V445 Lyr and RR Lyr itself with their values of 
0.0224 (= 16.6 min) and 0.0138 (= 1 1.3 min), respectively. There is 
no clear indication for any connections between the strengths of the 
two type of modulations. 



As was shown by ISzeidl & Jurcsik I 1 120091) and iBenko et al. I 

C20Q9), phase modulation always causes multiplet structures of 
higher order than triplets in the Fourier spectrum around the main 
frequency and its harmonics. These multiplet peaks were most 
clearly detected for V808 Cyg, RRLyr and V360 Lyr, the Blazhko 
stars with the strongest phase modulation in our sample. 

The precise and almost continuous observation of 14 Blazhko 
RR Lyrae stars measured by Kepler convincingly demonstrates that 
in all Blazhko stars both amplitude and phase variation are present. 



3.3 Additional frequencies 

Besides the modulation components that occur in multiplet struc- 
tures around the main frequency and its harmonics, we found addi- 
tional frequencies. In the cases of V808 Cyg, V355 Lyr and RRLyr 
these frequencies are located around /o/2, 3/o/2, 5/n/2, . . . , 
where /o denotes the main pulsation frequency. This very inter- 
es ting period doub l ing eff ect has already been discussed briefly 
inlKolenberg et al. 1 1 120 lOah for RR Lyr itself, and a separate paper 
dSzabo et al. 1120101) is dedicated to it. Here we just remark that the 
presence of these frequencies in the spectra seems to be variable in 
time and connected to particular Blazhko phases. The phenomenon 
can be described in a purely radial framework of pulsation theory. 

Similarly, time-dependent phenomena might also be impor- 
tant for the four further stars V354Lyr, V2178 Cyg, V360Lyr and 
V445 Lyr, where we discovered additional frequencies with small 
amplitudes. The light curves of V354 Lyr and V360 Lyr are plotted 
in the top left panels of Fig.[5] The last column of Table|2]indicates 
the possible identification of the additional frequencies. 

In the Fourier spectrum of V354 Lyr (KIC 6183128) the high- 
est (cr s = 34.4) additional peak is at f 2 = 3.0369 ± 0.0002 d" 1 
(see insert in Fig.[5}. Its ratio to the fundamental frequency (/n = 
1.78037 ± 0.00004 d -1 ) is 0.586, which is close to the canoni- 
cal ratio of the fundamental and second radial overtone modes. 
Furthermore, many linear combinations in the form kfo ± /2, 
k = 0, 1, 2, . . . are present in the spectrum. As the period dou- 
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Figure 6. Linear pulsation models of V350 Lyr and V354 Lyr. The observed 
frequencies were fitted with the theoretical fundamental and second over- 
tone modes. The lines represent the mass-luminosity values of the models 
which exactly fit the observations. Thin (red) lines: V354 Lyr periods with 
mefallicities (left to right): Z=0.0001, 0.0003, 0.001; thick (black) lines: 
V350 Lyr periods with metallicities (left to right) Z=0.0001 ; 0.0003; dashed 
(blue) line: V350Lyr periods with Z=0.0001 and the period ratio decreased 
by 0.001 to show the effect of possible non-linear period shifts. The thicker 
parts of the curves indicates the parameter range where the second overtone 
is linearly unstable. 



bling effect has a time-dependent nature, we have checked whether 
the mode connected to f 2 is also temporally excited or not. We 
tested this possibility both with the analytical function method and 
the amplitude variation tool of PERIOD04. These two independent 
methods yielded very similar results and showed the amplitude of 
f 2 changing over the Blazhko cycle (see the plot below the light 
curve of V354 Lyr in Fig.[5}. This means that V354 Lyr is a double 
mode pulsator with Blazhko effect, but not in a traditional sense. 

The observed frequencies /o and f 2 can be easily matched 
by the fundamental and second overtone modes in lin ear pulsation 
mode l s (Fig.[6]>. We used the Florida-Budapest code l lYecko et al.l 
Il998l ; iKollath et al. ||2002|) to fit the observed frequencies by the 
theoretical fundamental and second overtones. The model frequen- 
cies depend on four parameters, the mass, luminosity, effective tem- 
perature and chemical composition. Fixing two frequencies reduces 
the unknown parameters by two. Then for given chemical com- 
position and effective temperature the mass and luminosity of the 
matching star can be determined by numerical model calculations. 
The resulting mass-luminosity pairs are displayed on the Fig.[6] In 
the figure the linear instability region of the second overtone is also 
displayed. We have to note, however, that the width of these regions 
strongly depends on turbulent/convective parameters, and in addi- 
tion, the linear instability alone is not a sufficient condition for dou- 
ble mode pulsation. The coexistence of the fundamental and second 
overtone modes in non-linear pulsation is a new enigma that further 
theoretical studies should follow. 

The frequencies /' (2.0810 ± 0.0002 d" 1 ; a s = 28.1) /" 
(2.4407 ± 0.0003 d _1 ;cr s = 11.6) and /'" (2.6513 ± 0.0003 d" 1 ; 
<t s = 10.8, insert in Fig.[5jl are more difficult to identify. Their lin- 
ear combinations with /o and its harmonics also appear in the spec- 
trum and ratios are: f /f = 0.855, f /f" = 0.729, f /f" = 
0.672. The last two ratios permit an explanation as a first overtone 
(/" = /1) and a period-doubled (/'" = 3/o/2) frequency, re- 
spectively. However, /' does not seem to be fit this picture. This 
frequency might be connected to a non-radial mode. 

The star V2178 Cyg (KIC 3864443) shows a similar, but sim- 



pler and a bit less significant, additional peaks than V354Lyr. In 
this case the highest peak can be found again at the second over- 
tone frequency f 2 (3.5089 ± 0.0002 d _1 ; cr s = 15.8). One further 
peak at /' = 3.0585 ± 0.0003 d -1 (cr s = 9.8) was detected that 
might be a sign of a marginal period doubling effect (3/o/2). 

The other star where we found highly significant additional 
frequencies is V360Lyr (KIC 9697825). The dominating extra 
peak in its spectrum is located at /1 = 2.4875 ± 0.0001 d~ with 
<7 S = 51.9 (bottom right in Fig.[5}. The ratio of this frequency to 
the fundamental mode frequency (f = 1.79344±0.00003 d" 1 ) is 
0.721 . This value is a bit smaller than the generally accepted ratio of 
fundamental to first radial overtone frequencies (/1//0 = 0.745), 
but some model ca lculations with higher metallicity (see Fig. 8 in 
IChadid et al. 1120101) allow double mode pulsation with this ratio as 
well. We also detect a low-amplitude, but significant (<r s = 18.9), 
peak at the frequency of /' = 2.6395 ± 0.0002 d" 1 . Its ratio 
/o//' = 0.679 is almost the same as was found for V1127 Aql, 
a Blazhko RR Lyrae star observ ed with C0R0T and e xplained there 
by a non-radial pulsation mode ( IChadid et al. II 20 1 Of) . 

V445 Lyr (KIC 6 1 86029) has a very complex structure of fre- 
quencies at the low-amplitude level (see right panel in Fig.[2jl. Here, 
the three types of frequency patterns (connected to the first and 
second overtone modes and period doubling effect) show simi- 
larly significant peaks. In decreasing order of amplitude: 3/o/2 = 
2.9231 ± 0.0002 d^ 1 (a s = 24.5), /1 = 2.7725 ± 0.0002 d" 1 
(a s = 21.1) and f 2 = 3.331549 ± 0.0002 d" 1 (cr s = 13.5). 

If the global physical parameters are varying o ver the Blazhko 
cycle as demonstrated by ljurcsik et al. I j2009al lbl) excitation of a 
radial (or maybe a non-radial) mode temporally, (when the physical 
conditions of excitation are realized) would be a natural explanation 
for these patterns. 

The frequency spectra of non-Blazhko stars were checked 
to see whether any of them possesses additional frequencies. 
This search resulted in the discovery the double mode nature 
of V350Lyr (KIC 9508655). Its Fourier spectrum (see Fig.Q 
contains a significant (<x s = 11.5) peak at f 2 = 2.8402 ± 
0.0003 d" 1 with the amplitude of A(f 2 ) = 0.001 mag. The ratio 
of this frequency to the fundamental mode frequency (1.682814 ± 
0.00003 d _1 ) is 0.592. That is, V350Lyr is the first example of a 
non-Blazhko double mode RR Lyrae star, where the fundamental 
and second overtone modes are excited. The extremely high ampli- 
tude ratio, A(fo)/A(f 2 ) = 316, points out why we have not been 
able to find such stars from the ground. 

The residual peaks around the main frequency and its harmon- 
ics raise the possibility of the presence of the Blazhko effect close 
to the detection limit. However, neither the A\ (i) function nor the 
amplitude of the overtone frequency show clear time variation. In 
the future, calibrated Kepler data gathered over a longer time span 
may clarify the nature of this interesting object. 

In the case of V350 Lyr, the match of the observed frequen- 
cies to the theoretical model is a tougher problem than the case of 
V354 Lyr. High mass or luminosity is required to fit the empirical 
data even at low metallicity (Fig.[6). However, a small shift in the 
period ratio results in a better agreement with the canonical mass 
and luminosity values. 0.001 difference in the observed period ra- 
tio and the linear model values can be accounted to observational 
errors and non-linear effects of stellar pulsations. 
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Figure 7. Top: Phase diagram of V350Lyr (KIC 9508655). The Kepler 
light curve is folded by the main period of Po = 0.59424 d. Bottom: 
Fourier amplitude spectrum of V350Lyr around the highest peaks after the 
data were prewhitened with the main frequency /o and its harmonics. Due 
to small instrumental trends some residuals are seen at the location of re- 
moved harmonics (shown by green arrows). 



4 SUMMARY 

In this paper we have outlined some results obtained for RR Lyrae 
stars based on the first 138-d long data sets of the Kepler Mission: 

• We have determined the main pulsation periods and ampli- 
tudes for all stars in this sample. These parameters were previously 
unknown or wrong for 9 stars. 

• We have found 14 certain Blazhko stars among the 29 stars of 
the observed sample (48%), and we could not find any deviations 
from monoperiodicity for 15 stars (52%). Statistical distribution of 
the measured modulation amplitudes proved to be highly depen- 
dent on the used size of bins. Our sample supports the uniform 
distribution. 

• The possibility of multiple modulations of an RR Lyrae star 
is well known, such as the 4-year cycle of RR Lyr itself. Further 
secondary modulation cycles were also reported in the literature. 
We found a long term variation of the Blazhko cycle in V445 Lyr. 
The forthcoming long time base of the Kepler data will allow us to 
investigate this phenomenon in detail. 

• The Kepler data made it possible to find the Blazhko mod- 
ulation of KIC 1125706 with an amplitude as small as 0.03 mag. 
This is by far the smallest modulation amplitude ever detected for 
a Blazhko star. The same is true for phase modulations: we found 
small but clear phase variations for stars V2178Cyg, V783Cyg, 



V349Cyg, and KIC 1125706. For these cases the 8P < 1.5 min 
period variation again have the smallest known values. 

• The sensitivity both in amplitude and phase made it possible to 
detect phase variation for all Blazhko stars. Although our sample is 
small, it is notable that all of our Blazhko stars are modulated both 
in their amplitude and phase simultaneously. The relative strength 
of the two types of modulation varies from star to star, but always 
has a common period. Therefore, a plausible explanation for the 
Blazhko effect must account for both amplitude and phase varia- 
tions. 

• We found additional frequencies, beyond the main frequency, 
its harmonics and expected modulation components, in the Fourier 
spectra of 7 Blazhko type stars. These additional frequencies con- 
centrate around /o/2, 3/n/2, ... for three stars, while they appear 
around the first and second overtone frequency and its linear com- 
bination with the harmonics of the main frequency in the the case 
of V354Lyr, V2178Cygand V360 Lyr, respectively. A special case 
is V445 Lyr where all of the above-mentioned types of frequencies 
are present. 

If the basic physical parameters, such as mean radius, luminos- 
ity, effective temperature, are varying over the Blazhko cycle, a ra- 
dial (or maybe non-radial) mode could be excited temporally. This 
would be a natural explanation for these transients. 

• As a by-product of our frequency search for additional fre- 
quencies, we may have found the first double mode RR Lyrae star, 
V350 Lyr, which pulsates in its fundamental and second overtone 
mode simultaneously. 

Comprehensive and more detailed study of the increasing Kepler 
RR Lyrae data sets are in progress and will be discussed in the near 
future. 
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